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Abstract —In this paper, we propose a new coding scheme 
for symmetric Gaussian interference channels with feedback 
based on the ideas of time-varying coding schemes. The 
proposed scheme improves the Suh-Tse and Kramer inner 
bounds of the channel capacity for the cases of weak and not 
very strong interference. This improvement is more significant 
when the signal-to-noise ratio (SNR) is not very high. It is 
shown theoretically and numerically that our coding scheme 
can outperform the Kramer code. In addition, the generalized 
degrees-of-freedom of our proposed coding scheme is equal 
to the Suh-Tse scheme in the strong interference case. The 
numerical results show that our coding scheme can attain 
better performance than the Suh-Tse coding scheme for 
all channel parameters. Furthermore, the simplicity of the 
encoding/decoding algorithms is another strong point of our 
proposed coding scheme compared with the Suh-Tse coding 
scheme. More importantly, our results show that an optimal 
coding scheme for the symmetric Gaussian interference chan¬ 
nels with feedback can be achieved by using only marginal 
posterior distributions under a better cooperation strategy 
between transmitters. 

Index Terms —Gaussian Interference Channel, Feedback, 
Posterior Matching, Iterated Function Systems. 

I. Introduction 

The capacity of the interference channel with feedback 
has been still unknown for many decades although there 
were some progresses toward solving this problem. Kramer 
developed a feedback strategy and derived an outer bound 
of the Gaussian channel 16], 19]. However, the gap between 
the outer and the inner bounds becomes large unboundedly 
as signal to noise ratio (SNR) and interference to noise 
ratio (INR) increase. Furthermore, Suh and Tse 11], 12] 
characterized the capacity region within 2 bits/s/Hz and 
the symmetric capacity within 1 bit/s/Hz for the two-user 
Gaussian interference channel with feedback. They also 
indicated that feedback provides multiplicative gain at high 
SNR. However, their coding scheme does not work well 
when the SNR is close to the INR. Its symmetric coding 
rate becomes less than the Kramer code when this condition 
happens. In addition, it also has lower performance than the 
Kramer code when a = log INR/ log SNR is not very large 
and the SNR is low (cf. Figs. 1-2 of this paper or Fig. 14 
in 12]). Later, the Suh-Tse coding scheme was extended to 
M-user Gaussian interference channels with feedback for 
M > 3 17]. 

In this paper, we propose a coding scheme which 
achieves better performance than the Suh-Tse code when 
a = log INR/log SNR is not very large (See Figs.l- 
2 of this paper). In addition, our code can attain better 
symmetric rate than the Kramer code for all channel 
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Figure 1. Symmetric Rate Comparision at High SNR 
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Figure 2. Symmetric Rate Comparision at Low SNR 

parameters, and therefore it overcomes all the weak-points 
of the Suh-Tse coding scheme and improves the Suh- 
Tse and Kramer inner bounds. For the strong interference 
case, our code can achieve the same generalized degrees- 
of-freedom as the Suh-Tse coding scheme. Furthermore, 
numerical results show that our coding scheme indeed 
has better/equal performance than/to the Suh-Tse coding 
scheme for all channel parameters. 

H. Channel Model 

We consider the Gaussian interference channel shown in 
Fig. 3, which has two senders and two receivers. Sender 1 
sends a source of information message points ©i, which 






































Figure 3. Gaussian Interference Channel with Feedback 

is uniformly distributed in (0,1), to receiver 1. Sender 2 
sends a source of information points © 2 , which is also 
uniformly distributed in (0,1), to receiver 2. Assuming that 
©1 is independent of ©2 and each channel interferes with 
the other. Specially, we assume that 

Fi =Xi+aA2 + Zi, 


The rate pair is achieved within input power constraints 
Pi, P 2 if the following is satisfied: 

limsup - < Pm^'rn G {1,2}. (2) 

The symmetric capacity is defined by Csym •= supjP : 
(P, P) is achievable}. 

An optimal fixed rate decoding rule for the two-user 
Gaussian interference channel with feedback for rate pair 
(Pi,P 2 ) is the one that decodes a pair of fixed length 
intervals {(Ji,J 2 ) : \Jm\ = for m G {1,2}}, 

which maximizes posteriori probabilities, i.e., 

AM(y(n.m)) = argmax 

It is easy to see that the optimal fixed rate decoding rule 
for the Gaussian interference channel with feedback is the 
traditional MAP, MMSE decoding rule. 

An optimal variable rate decoding rule with target error 
probabilities (n) = 6^^ is the one that decodes a pair 
of minimal-length intervals (Ji, J 2 ) such that accumulated 
marginal posteriori probabilities exceeds corresponding tar¬ 
gets, i.e.. 


I2 — ^2 P aXi + ^ 2 , 


where ^ N'{0,al) and Z 2 ^ A/'(0,cr2) are Gaussian 
noise random variables, and the input power constraints 
are Pi and P 2 , respectively. In this paper, we consider the 
symmetric interference channel such that Pi = P 2 = P 
and (Ji = (72 = 1. Hence, the signal-to-noise ratio and 
interference-to-noise ratio can capture channel gains SNR 
= P and INR = a^P. We also assume that output symbols 
are casually feedbacked to the corresponding sender and 
the transmitted symbol at time n can depend on both 
the message ©^ and the previous channel output sequence 
Yin-i,m) for m e {1,2}. 

A transmission scheme for the two-user Gaussian inter¬ 
ference channel with feedback is sequences of measurable 
functions : (0,1) x ^ {1^2} so 

that the input to the channel generated by the transmitter 
is given by 




A decoding rule for the two-user Gaussian interference 
channel with feedback are sequences of measurable map¬ 
pings G {1, 2} where S is the set 

of all open intervals in (0,1) and refers to 

the decoded interval at receiver m. The error probabilities 
at time n associated with a transmission scheme and a 
decoding rule, is defined as 


e {1,2}, 

and the corresponding rate pair {Rn\Rn^) at time n is 
defined by 

ijW ;= -1 log |aM 

We say that a transmission scheme together with a de¬ 
coding rule achieves a rate pair (Pi,P 2 ) over a Gaussian 
interference channel if for m G {1,2} we have 

lim P = 0, lim = 0. (1) 

n^oo \ J n^oo 


AM(yK-))= min |J„|. 

Both decoding rules use the marginal posterior distribution 
of the message point P©^|yrr which can be calculated 
online at the transmitters and the receivers. Refer [6] for 
more details. 

Lemma 1: The achievability in the definition (1) and (2) 
implies the achievability in the standard framework. 

Proof: See the detailed proof in [3], [4], and [5]. ■ 

Notations: The cumulative distribution function (c.d.f.) 
of a random variable X is given by Fx{x) = 
IPx((—oc,x]), and their inverse c.d.f. is defined as 
(t) := mf{x : Fx{x) > t}. The composition 

function is defined by (/ o g){x) = f{g{x)), and the 
sign function is defined as sgn(x) := 1 if x > 0 and 
sgn(x) := —1 if X < 0. We also use (x)+ := max(x, 0) 
and log^(x) := max(logx,0). 

III. A Coding Scheme for Gaussian 
Interference Channels with Feedback 

In this section, we propose a time-varying coding scheme 
for the symmetric Gaussian interference channel with feed¬ 
back as following: 

A. Encoding 

m Step 1: Transmitter m sends x[^^ = P^^(©^) 
where m G {1,2} and X ^ A/’(0,Pi) for some 
Pi > 0. We also set 

Pi .=--=0. 

• Step n + 1 for n > 1: 

Both transmitters estimate 

Pn+l = {pn ~ 2^nSgn(/9^)(|/9^| + |a|) 

Pn 

+6^sgn(/)n)[|/)n|(l + |ap) -f 2|a|]} . 



















Transmitter 1 sends X^^^sgn(/)^+i) where 

:= - &nsgn(p„)y;(i)). 

Pn 

Transmitter 2 sends X^^^sgn(a) where 

:= T(x( 2) - &„sgn(a)F„(2)), 

Pn 

Receiver 1 receives 

= X«iSgn(p„+i) + \a\xl^P + zp,. 

Receiver 2 receives 

= sgn{a)xl^P + a sgn(p„+i)X«i + 

Both receivers feedback their received signals to the 
corresponding transmitters. 

Here, {/3n > 0,6^) should be chosen to satisfy the 
following constraints: 

1 ^ 

limsup - V E[Xpp < P, Vm e {1,2}. (3) 

iV^oo 

B. Decoding 

• At each time slot n, receiver m G {1,2} selects a 
fixed interval = {sm^tm) C M as the decoded 
interval with respect to X^\. 

• Then, set the decoded interval = 

as the decoded interval 

with respect to x[^\ where 
T^^\s) := o o • • • o w^\s)^ Vs G M 

and 

wP ■■= /3nX + bnSgn{pn)Yp\ 
wP := PnX + bnSga{a)Yp\ 

• Receiver m sets the decoded interval for the message 
0^ as follows: 

('yK-)) = FxMP^)- 

We call this coding strategy the Gaussian interference 
time-varying feedback coding strategy, which is an optimal 
variable rate decoding rule with doubly exponential decay 
of targeted error probabilities (see the proof of the Lemma 
2 in this paper). 

IV. A New Achievable Rate Region 

Lemma 2: Under the condition that 0 < 

lim sup^^QQ < 1, the time-varying coding scheme 

for the symmetric Gaussian interference channel with 
feedback achieves the following symmetric rate: 

Rsym = — hm sup log f3n (bits/chanuel use). 

n^oo 

Proof: We provide a sketch of the proof. The detailed 
one can be found in papers [3], [4], and [5]. 

• Define {3 := limsup^^o^ It is easy to see that 

^sym — ^ 

m For any R < Rsym, we can find an e > 0 such that 
R < \og{P + e)“V 

• Choose an Vg G N such that sup^> < P -\- e. 


• Using the Law of Iterated Expectations (Fubini’s 
Theorem), we can show that 

V{Rp'> <R)< + 

• We can also show that 

By choosing ~ 

^ l^ 2 '^(Rsym-R)^^ and tm = —Sm = \j[^"^\/2, the 

aforementioned condition (1) is satisfied. Here, sym¬ 
bols 0{x) and o{x) are Landau symbols. Besides, the 
choice of sequences {bn,(3n) following the rule (3) 
leads to the fact that the input power constraints are 
also satisfied. 


Theorem 1: The non-degraded symmetric Gaussian 
interference channel (a 7 ^ 0 ) can achieve the following 
symmetric rate: 


^sym(bits/channel use) = - 


max 


X log 


2 pe[0,po], be{bl,b2} 

p 


P + 62[1 + P + |app -f 2|a|Pp] - 26P[1 -f |a|p] 


where 


/a2p2 + P - v'P[2a2p2 + Pl 

0 < Po y < 1 , 


and 


K 


2Pp+\a\P+\a\Pp^ 


^ 2|a|P + 2Pp + 2|a|2Pp + p + 2|a|Pp2 

^ yp2|a|2p4 _ 2p2(|a|2p2 + P) + |a|2p2 
^ 2|a|P + 2Pp + 2|a|2Pp + p + 2|a|Pp2 ’ 

Proof: From our transmission strategy, we have 

xp, = Xxp - bnSgn{p„)YyP (4) 

Pn 

Xpl, = f(XP - 6„sgn(a)yi2)). ( 5 ) 

Pn 


Denote 


P„ = E[XP]^ = P[X(2)]2,p„ = ^ p " ^ 

-Rn 

(We can show by induction that E[xPY = ElxP]"^ for 
all n). Therefore, it is easy to see that 

E[XPYW]=PM + \a\), 

E[xPyJiP = P„sgn(a)sgn(p„)(|p„| + |a|), 

ElYp^Yp^] = P„sgn(a) [|p„|(l + \a\^) + 2\a\] . 

Note that 

E[xp,xp,] = T t^E[XPxP] 

-6„sgn(p„)P[x7)Fii)] - 6„sgn(a)P[x7)F„(2)] 

+^nSgn(p„)sgn(a)P[r7)y7^]} • 



Finally, we obtain 


Pn-\-lPn-\-l — Pn^E^iPn) 02 ‘^^n{\Pn \ “1“ |^|) 

Pn 

+^n[|Pn|(l + hl^) + 2|a|]} • (6) 

Similarly, observe that 

^[xWrii)] = P„sgn(p„)[l + |a||p„|], 
^[Xi2)ri2)] = P„sgn(a)[l + |a||p„|], 

= l + Pn + a^Pn + 2\a\\pn\Pn, 


Note that since Lemma 2 holds only for 0 < < 1, the 

superscript + is necessary to deal with general cases. 

To complete the proof, we show a procedure to force 
\Pn\ = p{pn = (-l)"/o) for any p e [0,po], K = b, 
Pn = P, Pn = P for any n> 2. Indeed, from (6) and (7), 
we firstly force P 2 = P^ p 2 = p, and P 2 = P setting 

Pp=^^{2\a\{bi-b^)}, (11) 

P = 4{-Pi - 261 P 1 + blil + Pi + a^Pi)}. (12) 

Pi 


P[F„(2)]2 = 1 + P„ + + 2|a||p„|P„. 

From the relations (4) and (5) we have 

Pfi+l = {Pfi — 2Pn^n[l + l<^l|Pn|] 

Pn 

+6^ [1 + Pn + + 2|a| |/9n|Pn] } • (7) 

From ( 6 ) and (7), if we can force P^ P^bn Pn ^ P 
and |pn| ^ P ^ [O 5 1]’ we obtain the following equations 
with three unknowns {b, p, P)\ 

P=P[P- 26P(1 + \a\p) + 6"(1 + P + a"P + 2\a\pP)] , 

(8a) 

—/?= — [/? — 2b{p + |a|) + (/3(1 + |a|^) + 2|a|)] . (8b) 

By eliminating p and considering p as a running variable, 
we have the following quadratic equation in b for each 
fixed choice of p\ 

6^[2|(2|P -|-2Pp-j-2|(2|^Pp-|-P-|-2|(2|Pp^] 

- 26[2Pp + |a|P + P\a\p‘^] + 2Pp = 0. (9) 


After some simple derivations, the discriminant of this 
quadratic equation is given by 

A = p2|a|2/-2p2(|a|2p2^p)^|^|2p2 _ 

Since /(O) = japP^ > 0 and /(I) = —2P < 0, there 
exists the minimum value po ^ (0? 1) such that /(po) = 0- 
Specifically, the value of po is given by 


/ |a|2p2 + p - JP[2a‘^P‘^~Ty] 

0 < /’o - y < 1. 

On the other hand, since the derivative of f{p) satisfies 
/'(p)=4p2|a|2(p3-p)-4Pp<0, 


for all p G [0,1), we have A = f{p) > 0 for all p G [0, po]- 
For all these values of p, we can easily show that (9) can 
have two positive solutions b ^, b^ as the theorem statement. 
As a result, (8a) and (8b) have at least one solution (6, P) 
for each fixed p G [0,po]- Therefore, if we force all the 
sequences {bn, Pn,PnApn\) to converges to {b,P,P,p), 
the symmetric Gaussian interference channel with feedback 
can achieve the following rate by the Lemma 2: 


Psym — 


— log( min P) 

p^[0,po] 


1 

- max X 
2 pe[o,po],be{bi,b^} 


This procedure is feasible because (11) and (12) have at 
least one solution which is a triplet (6i,Pi > 0, Pi > 0) 
for each p G [0,po]- Indeed, 

• For p = 0, we can choose bi = 0, Pi = P, pi = 1. 

• For p 7 ^ 0, from (11) and (12) we have the following 
quadratic equation in bi 

6^[(1 + Pi +a^Pi)p — 2|a|Pi] — 2(p— |a|)Pi6i +Pip = 0. 

(13) 

The discriminant of this quadratic equation can easily be 
shown to be equal to P(p) = a^(l — p^)Pf — Pip^. For the 
case p / |a|, we can choose Pi such that this discriminant 
is equal to zero by setting Pi = p^ j{oP(\ — p^)). By this 
choice of Pi, we obtain 

{p-\a\)Pi 

^ (l + Pi+a2Pi)p-2|a|Pi' 

In order for (11) and (12) to have solution pi, we need to 
show that the above choices of Pi, 6i satisfy b^ — bi > 0. 
Clearly for p < |a|, this requirement is satisfied by noting 
that 6i < 0 since p — |a| < 0 and 

^ 2|a|p2 ^ 2|a|Pi 
^ a^+p^ 1 + Pi + a^Pi ■ 

For p> \a\ (|a| < 1, of course), observe that 

p^ ^ p 

^ a2(l-p2) > \a\-a^p- 

It follows that (p—|a|)Pi > (1 + Pi+a^Pi)p —2|a|Pi >0 
or 6i > 1. This also means that bl-bi> 0. For the case 
p = |a| any choice of Pi > 1/(1 — a^) works since we 
have E{\a\) > 0. Besides, the sum of two solutions of the 
quadratic equation (13) in bi is equal to zero, and their 
product is not equal to zero (by using Vieta’s formula). 
Hence, we must find at least one 6i < 0 or 6^ — 6i > 0. 

Finally, we only need to set pn = p,Pn = P, bn = b 
which is a solution of (8a) and (8b) for each choice of 
p G [0,po] and for all n > 3. Here b should be chosen 
to minimize p for each choice of p G [0,po] in order 
to maximize the achievable symmetric rate of our coding 
scheme. Last but not least, we can show that our code has 
better performance than the Kramer code [6], and therefore 
the superscript + can be got rid of from the achievable 
symmetric rate formula. ■ 

Remark 1: For the degraded Gaussian interference chan¬ 
nel with feedback (a = 0), (8a) and (8b) become 


log 


P + 62[1 + P + \app + 2|a|Pp] - 26P[1 -h |a|p] 


-P= ^P{b - l)^ P = T[P(6 - 1)2 + b\ (14) 



0 and the achievable 


From (14), we must have p = 
symmetric rate becomes 


= llog(l + P). 


Rsym = o “-a? log 


beR 


P( 6 - 1)2 + 62 


positive number ^ < otj2 although this p^ may not be 
optimal. Indeed, for P sufficiently large, we have p* ~ 1 
and 

p* < - p-i + y^p-io>-i) -p-io^-i)/2 = y/i-p-i. 


This result coincides with the well-known capacity of this 
channel with no interference. 

Corollary 1: The proposed time-varying code outper¬ 
forms the Kramer code for all channel parameters. 

Proof: A variant of Kramer code is constructed by 
choosing triplet which is a solution of ( 8 a) and 

( 8 b), as follows: 


P (1 + \a\p) I ffl 2 p(i_p 2 )+~ 

P(l+a2 + 2|a|p) + l’^ y P(l+a2 + 2|a|p) + l’ 

and p is the unique solution in ( 0 , 1 ) of the next equation: 

2\a\^P‘^p^ + a^Pp^ - 4\a\P{a^P -f l)p^ 

- {2a^P + P + 2)p + 2|a|P(a^P + 1) = 0. (15) 


(See also in [2], [ 6 ], [9]). Therefore, it is inferior to the 
proposed code in this paper. 

Remark 2: The choice of b and (3 for the Kramer code is 
to maximize the achievable rate (or minimize /3) by using 
only ( 8 a) for each fixed value of p. In order to satisfy ( 8 b), 
the choice of 6 , p in the above proof is applicable only 
to the fixed value of p which is the unique solution in 
(0,1) of (14). For other values of p, from ( 8 a) and ( 8 b), 
we see that b and [3 are given by two different functions 
of p. Enlarging the set of possible choices of p (then 6 , /3) 
increases the achievable symmetric rate in this paper. ■ 

Corrolary 2: For a = log/A^P/log > 1, the 

generalized degrees of freedom of the proposed coding 
scheme is given by 


p(^\ ._ 1 ™ RsymjSNR, I NR) _ a 

SNRJNR^oo logPA^P 2' 

(Here, the unit of Rsym is bits/s/Hz). As a consequence, the 
proposed coding scheme has the same generalized degrees 
of freedom as the Suh-Tse coding scheme [1], [2] and also 
provides the multiplicative gain at high SNR. Note that the 
Kramer code achieves only d{a) = (1 + a)/4 for a > 1 
(cf. (48) [2]). 

Proof: Since a = log/A^P/logPA^P > 1, we have 
|apP = P^, or \a\ = p(«-i )/2 por P sufficiently large 
and p ~ 1 , observe that 


p(a+l)/2(l + p(a+l)/2(^ _ ^2) 




6 * - ^ 

’ 2 P«p 2 P^p 

By choosing b = b^ ^ we obtain 

Rsymip) ~ log (p ^ ^2pa _ 2bP{l + P(a-l)/2p) 


= - log (l - p - (bits/s/Hz). (16) 

From Theorem 1, we also have 


P = 1 + p-« - v' 2 P-« + P- 2 « > 1 - V3P-°‘/^, 


and the p*, that maximizes the achievable rate, must be 
in the interval [0,po]. This condition is satisfied by setting 
p^ = Vl — P-'^ -f p-(«-i) — p-(a-i )/2 fQj. arbitrary 


Hence pj<l — P ^<1 — ^/3P < pg. On the other 

hand, since 1 — pi — 2 p*P^^“^^/^ = P“^, we obtain 


d{a) > lim 

SNR,INR^oo logPA^P 


Since 7 can take any arbitrary value less than a/2, we 
have d{a) > a/2. From the result of [2], it is known that 
d{a) < a/2. Hence, we must have d{a) = a/2. ■ 


V. Numerical Evaluation 

In order to evaluate Rsym in Theorem 1 numerically, 
we determined the optimal p by increasing it from zero by 
incremental step 10“^. The numerical results are shown 
in Figs. 1-2 in Section I. We note that our proposed code 
can achieve better/equal symmetric rate than/to the Suh- 
Tse and Kramer codes for all channel parameter (a,P). 
This improvement is more significant when SNR is not 
too high. 


VI. Conclusion 

The inner bound of the capacity is improved compared 
with the Suh-Tse and Kramer inner bounds by analyzing 
the performance of our proposed code as an optimized 
form of the Kramer code. Our result also shows that 
an optimal coding scheme for the Gaussian interference 
channel with feedback can be achieved by using marginal 
posterior distributions. 
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